Effective temperature for hopping transport in a Gaussian DOS 
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For hopping transport in disordered materials, the mobility of charge carriers is strongly de- 
pendent on temperature and the electric field. Our numerical study shows that both the energy 
distribution and the mobility of charge carriers in systems with a Gaussian density of states, such as 
organic disordered semiconductors, can be described by a single parameter - effective temperature, 
dependent on the magnitude of the electric field. Furthermore, this effective temperature does not 
depend on the concentration of charge carriers, while the mobility does depend on the charge carrier 
concentration. The concept of the effective temperature is shown to be valid for systems with and 
without space-energy correlations in the distribution of localized states. 

PACS numbers: 72.20.Ht, 72.20.Ee, 72.80.Le, 72.80.Ng 
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I. INTRODUCTION 

For many years much attention of researchers has been 
devoted to the effect of high electric fields on the hopping 
mobility of charge carriers in organic disordered materi- 
als, such as conjugated and molecularly doped polymers 
and semiconducting organic glasses Ja^i^i^ Two peculiari- 
ties in the dependence of the carrier mobility on the elec- 
tric field are usually discussed in the literature. One of 
them is the apparent decrease of the mobility with rising 
field at relatively low fields and high temperatures re- 
ported in several timc-of-fiight studies. ^'^iSi"^ The other 
one is the very strong non-linear increase of the car- 
rier mobility with electric field at low temperatures and 
high fields.— >^i^ While the former peculiarity has been 
attributed^i^ to misinterpretation of experimental data, 
the latter one has been confirmed in numerous experi- 
mental studies and it is currently in the focus of intensive 
theoretical research. Before we turn to discussing this 
research, it is worth noting that very similar non-linear 
effects with respect to the applied electric field have been 
known since decades for transport phenomena in the inor- 
ganic noncrystalline materials, such as amorphous semi- 
conductors. Indeed, strong nonlinearities in amorphous 
semiconductors were observed for the field dependence of 
the dark conduct ivityfSii.° of the photoconductivityjii and 
of the charge carrier drift mobility^ d^d'^ at high electric 
fields. While the field-dependent hopping conductivity at 
low temperatures was always a challenge for a theoretical 
description, the theories for temperature dependence of 
the hopping conductivity at low electric fields were suc- 
cessfully developed for all transport regimes listed above 
(see for instance Ref. y and references therein). Further- 
more, it has been shown that the effect of a strong elec- 
tric field on transport coefficients in amorphous semi- 
conductors can be effectively described by replacing the 
temperature parameter in the formulas for the low-field 
temperature-dependent mobility and conductivity by an 
effective temperature Tcs(T, F) , dependent on the mag- 
nitude of the electric field F. Shklovskii^'' was the first 



who recognized that, for hopping conduction, a strong 
electric field plays a role similar to that of temperature. 
In the presence of the field, the number of sites avail- 
able for charge transport is essentially enhanced in the 
direction prescribed by the field. The distance between 
sites available for hopping transport shortens and hence 
electrons can move faster. -'^'^ The concept of the effective 
temperature has been studied in detail for systems with 
an exponential energy distribution of localized states usu- 
ally assumed for inorganic noncrystalline materials: 
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where N is the concentration of localized states and a 
is the energy scale of the distribution. By studying the 
steady-state energy distribution of electrons in numerical 
calculations and computer simulations^^"^^ and by com- 
puter simulations of the steady-state hopping conductiv- 
ity and the transient energy relaxation of carriersi^ the 
following result has been found. The whole set of trans- 
port coefficients can be represented as a function of a 
single parameter reff(r, F): 
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with (3 — 2 and values of 7 in the range 0.5 - 0.9 depend- 
ing on which transport coefficient is consideredj^i In this 
expression a is the localization length of charge carriers 
in the localized states, e is the elementary charge, and 
k is the Boltzmann constant. Herewith the problem of 
nonlinearities of transport coefficients with respect to the 
applied electric field for inorganic noncrystalline materi- 
als with the density of states (DOS) described by Eq. ([T|) 
has been solved. 

Let us now turn to the organic disordered materials. In 
such systems the density of states involved in the hopping 
transport of charge carriers is believed to be not purely 
exponential as in Eq. ([T]), but rather to be described by 
a Gaussian energy distribution>ii^i^i^ d^d^ 
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Numerous computer simulations have been devoted 
to studying the field nonlinearities of transport coeffi- 
cients in the hopping regime in such systems. Two mod- 
els of disordered organic materials were considered: the 
so-called Gaussian disorder model (GDM) suggested by 
Bassler et al.— and the so-called correlated disorder model 
(CDM) considered by Garstein and Conwell,— by Dun- 
lap et al-r^ and by Novikov et al^^^ In both models, 
the field dependence of carrier mobility has been studied 
by computer simulations. While analytical calculations 
have been carried out in order to justify the COM j^^i^^i^^ 
a consistent analytical theory for the field dependence of 
the hopping mobility in a Gaussian DOS is still missing. 

It is tempting to try to apply the concept of the effec- 
tive temperature to organic disordered systems with the 
Gaussian DOS described by Eq. ([3]), since this concept 
has proven to be very successful for inorganic systems 
with the exponential DOS described by Eq. ([1]). The 
idea that a strong electric field leads to heating of the 
charge carriers in organic materials has already been con- 
sidered in several theoretical studies. Recently Preezant 
and Tessler— performed such a study in the framework 
of an analytical approach that artificially decouples the 
energy-dependent factors in the hopping transport from 
the space-dependent factors. However, this approach has 
been shown2i2^ to be unsuitable for describing hopping 
transport processes in disordered materials. Moreover, 
such an approach leads to an effective temperature that 
essentially differs from the one considered in previous 
studies for inorganic materials^ii^ii^ii^iii where the ef- 
fective temperature arises from the interplay between 
spatial- and energy dependent factors in hopping pro- 
cesses. Li, Meller and Kosina3^ recently approached this 
problem by inserting the dependence of transition rates 
on the electric field and that of the percolation threshold 
into the percolation theory of Vissenberg and Matters^!. 
More recently, Limketkai et ali^ exploited the con- 
cept of the effective temperature in order to account 
for the strong field nonlinearity of the conductivity and 
carrier mobility observed in organic semiconductors i^ 
Limketkai et al. chose the expression for the effective tem- 
perature Tcs{T, F) in the form of Eq. ^ with j3 = 1 and 
7 = 0.5. It has been proven however^^ that such expres- 
sion for the effective temperature with /9 = 1 cannot be 
considered as suitable. Indeed, suppose the conductivity 
G is dependent on Tcff (T, F) solely. Then 
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In the Ohmic transport regime at i^ <C kT/ea, the con- 
ductivity G must be field independent, implying that 



dF 



The function described by Eq. ([2]) with /3 = 1 obviously 
does not fulfill this condition. However, any function of 
this kind with /3 > 1 is consistent with Eq. ([5|) along with 
the necessary condition Tcs ^ T at F = Q and Tcs oc F 
at T = 0. 

Therefore the challenging problem arises to find out 
whether the concept of the effective temperature is ap- 
plicable to systems with a Gaussian DOS and if yes, what 
is the expression for Tcff (T, F) . We try to answer these 
questions in the present study. For this purpose we fol- 
low the idea of Marianer and Shklovskii-'^^ to look at the 
energy distribution of charge carriers using the numeri- 
cal method of nonlinear balance equations suggested by 
Yu et ali^i^i^^ii^^ and apply it to electron transitions 
in a Gaussian DOS. Preliminary data for the energy 
distribution of charge carriers in the GDM at low car- 
rier concentrations confirmed the concept of the effective 
temperaturei'?* 

The paper is organized as follows. In Sec. [TT]the nu- 
merical method used in our study is described. In Sec. lIIII 
the results for the energy distribution function of charge 
carriers in the GDM are presented. It is shown that at 
finite temperatures and electric fields the energy distribu- 
tion of charge carriers in a Gaussian DOS can be well de- 
scribed by the Fermi-Dirac distribution function charac- 
terized by Toff(r, F) given by Eq. © with /3 = 1.54 ±0.2 
and 7 = 0.64 ± 0.2 thus confirming the applicability of 
the concept of the effective temperature to organic disor- 
dered solids. Furthermore, we show that the expression 
for Tcs{T, F) is stable against changes in the concentra- 
tion of charge carriers. Numerical results obtained for the 
mobility of charge carriers also confirm the concept of the 
effective temperature for the GDM. In Sec. II VI the corre- 
sponding results for the CDM are presented. They prove 
that the concept of the effective temperature is valid also 
for correlated systems. Concluding remarks are gathered 
in Sec. IVl 



II. NUMERICAL METHOD 

In order to check the validity of the concept of the 
effective temperature, we solved a system of non-linear 
balance equations with respect to the steady-state occu- 
pation probabilities pt of sites and checked, following the 
idea of Marianer and Shklovskii^ whether this distribu- 
tion can be fitted by a Fermi-Dirac distribution 
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with some single parameter Tctf{T,F). We also calcu- 
lated the carrier drift velocity along the field direction 
and the corresponding mobility at different temperatures 
and electric fields and checked whether the mobility can 
be described as a function of a single parameter Toff (r, F) 
combined from the temperature T and the field strength 
F. 



The balance equation for the occupation probability pi 
of a site i has the form^SiSiiHiS 

^p,ry(l -pj) = ^pjrj,(l -k), (7) 

where the rate of jumping from site i to site j is given by 
the Miller- Abrahams formula, 
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Here Ae^j is the difference between energies of states j 
and i, A/J^j is the distance between these states, and i^o 
is the attempt-to-escape frequency. While in the initial 
work of Marianer and Shklovskii^^ the limit of a low con- 
centration of charge carriers was considered and the bal- 
ance equations were linearized, we study a system with 
finite number of charge carriers n. In order to solve the 
system of the nonlinear balance equations we use the iter- 
ative numerical procedure suggested by Yu et al.^&'^'^i^ 
The balance equation ([7]) is rewritten in the form 
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where the right hand side does not contain pi . This ex- 
pression is used iteratively to find a solution to Eq. ([7]). 
In each iteration step, all pi are updated. If pj has al- 
ready been calculated in the current step, this value is 
used, otherwise the result from the previous step is used. 
This so called implicit iteration is necessary for obtaining 
convergence.— In the sums over j, only the most impor- 
tant transitions are needed. We have considered jumps 
shorter than a cut-off length i?cutoff, chosen so large that 
it does not affect the transport parameters. 

The procedure is repeated until the relative change of 
any one of the probabilities is smaller than 10~^°. We 
note that the solution seems to converge faster if the 
sites are placed on a lattice than at random. When the 
localization length is small, more iterations are needed, 
but the cut-off length can be reduced. 

If an electric field F is applied in z direction, the dif- 
ference in energy between sites j and i is given by 



Asij = Sj - £i - Fe{zj - Zi). 



(10) 



When the probability distribution of the charge car- 
riers over the sites is known, the drift velocity of the 
carriers along the field and their mobility are given by 

Vz^ Y P^^yd -Pj)(2;j - z,)/n, M = y • (H) 

We studied a system of M sites distributed randomly 
within a cube with a side length L. Each site has a 
random energy e, from a Gaussian distribution with the 



width (J. The density of states is given by Eq. ([3]). 
Periodic boundary conditions were applied in all direc- 
tions. The calculations were performed in dimensionless 
units, where the width of the DOS a, the site concentra- 
tion N — M/ L^ , the Boltzmann constant k, elementary 
charge e, and attempt-to-jump frequency vq are equal to 
unity. 

We used a system containing M = 8000 sites inside a 
cube with side length L = 20. The localization length 
was varied between 0.2 and 1, while the particle concen- 
tration was varied in the interval between c = 10^^ and 
c = 0.1. 

For each choice of T, F, concentration c, and local- 
ization length a, a number of realizations of the system 
were generated and the steady state occupation proba- 
bilities were determined. For each realization, the carrier 
mobility and effective temperature were calculated. The 
magnitudes of the effective temperaures obtained with 
different realizations of the system were equal to each 
other within the accuracy of one percent. 

The effective temperature was determined from the oc- 
cupation probabilities pi with a linear fit of ln(l/pi) as a 
function of e^. 
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This approach works well when the variation in occupa- 
tion probabilities for sites with similar energies is small. 
When the variation is larger, we found it better to pro- 
duce a histogram of p{e) and then fit to this histogram 
as above. The histogram approach seems more correct, 
since it uses the average occupation probability at a given 
energy, in contrast to the direct fit that uses the average 
of ln(l/p). All results for effective temperature presented 
below were obtained with the histogram method. The 
parameters P and 7 in Eq. ^ were determined by simul- 
taneous fitting of all data points to the surface given by 
this equation. 

In the next Section we present the results of calcula- 
tions for a system without any correlation between site 
energies and their spatial positions (GDM). In Section 
II VI we present the results of calculations for the CDM - a 
system, in which the energies of neighboring sites are cor- 
related. In order to introduce correlations, we generate 
an initial energy for each site with a Gaussian distribu- 
tion. We then calculate the real energy for each site by 
averaging the initial energies of all sites inside the corre- 
lation radius i?corr, in accord with the recepie suggested 
by Garstein and Conwell.^*^ This averaging procedure has 
two effects. It will introduce a correlation in energy be- 
tween spatially close sites and it will decrease the width of 
the energy distribution. In order to keep the width of the 
energy distribution independent of the correlation length, 
the site-energies were rescaled to provide the initial en- 
ergy width a. We discuss these effects of space-energy 
correlations in more detail in Section [IVI 

At very high electric fields, when the energy landscape 
is strongly sloped along the field direction and the energy 





FIG. 1: Occupation probability as a function of site energy, 
for three different charge carrier concentrations c. The tem- 
perature and field are T = 0.2 and F = 0.4, respectively. 
The dots show occupation probabilities of individual sites, 
while the symbols show a histogram, i.e. average occupation 
probability for sites in an energy interval. The curves show 
Fermi-Dirac functions fitted to the histogram. 
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FIG. 2: Calculated Tcfi as a function of charge carrier con- 
centration, for different values of T and F. 



differences between the successive states become negative 
as given by Eq. (fTO| . the carrier drift velocity saturates 
and becomes field- independent.^ We study the range of 
electric fields, vifhich are less than the field at which the 
velocity saturates. 



III. EFFECTIVE TEMPERATURE FOR 

SYSTEMS WITHOUT SPACE-ENERGY 

CORRELATIONS 

In Fig. [1] the numerical results are shown for the oc- 
cupation probability for each site, as a function of site 
energy, for three different charge concentrations. The in- 
set shows the same data on a logarithmic scale. 



FIG. 3: Tcff as a function of temperature T and aF. Con- 
centrations c — 10~ and f 0~ are shown. The surface is the 
best fit to Eq. Q, with j3 = f .54 and 7 = 0.64 



At low concentrations, the data can be fitted by 
a Boltzmann distribution with a parameter Toff (T, F) , 
while the data at high concentrations show a Fermi-Dirac 
shape, remarkably corresponding to the same parame- 
ter Toff(T, F). Furthermore, Fig. [5] shows that the value 
of Tofr(r, F) does not depend on the carrier concentra- 
tion c. These results are really remarkable because they 
mean that the concept of the effective temperature is 
valid for systems with a Gaussian density of states (so far 
it has been proven only for exponential DOS described 
by Eq. ([T])) and that the parameter Tcb{T, F) is universal 
with respect to the concentration of charge carriers in the 
system under study. 

In Fig. [3] we bring together all the calculation results 
obtained for various field strengths in the range < F < 
3 and for various temperatures in the range 0.1 < T < 0.5 
for concentration of carriers c = 10^^ and 10~^, in the 
form of a three-dimensional plot of Tcs{T, F), obtained 
by the best fit of the calculated energy distributions by 
a Fermi-Dirac function. Also the surface determined by 
Eq. ^ with parameters f3 = 1.54; 7 = 0.64 is shown in 
this figure. One can see an excellent agreement between 
the calculated results and Eq. ^ with these parameters. 

As described in the previous Section, we also stud- 
ied the mobility of charge carriers with respect to the 
applicability of the concept of the effective tempera- 
ture. In Fig. |4] the dependences of the carrier mobil- 
ity as a function of Tcf({T, F) determined from the fits 
of the energy distribution function are shown for dif- 
ferent concentrations c of charge carriers. These data 
were obtained for various field strengths in the range 
< i^ < 3 and for various temperatures in the range 
0.1 < T < 0.5. The data were averaged over three re- 
alizations. For the chosen parameters, it appears not 
important whether the mobilities or the inverse mobili- 
ties were averaged. While the magnitude of the mobil- 
ity appears sensitive to the concentration of carriers, in 
agreement with the results obtained previously by several 
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FIG. 4: The dependence of carrier mobility on the effective 
temperature, for different locahzation lengths a. In each pair 
of curves, the upper one is for charge carrier concentration 
c = lO"'^ and the lower one for c — 10~^. 



FIG. 5: Occupation probability in the correlated disorder 
case, with _Rcorr = 2. The temperature and field are T — 0.2 
and F = 0.4, respectively. The symbols are defined as in 

Fig.m 



research groups r^i21i2£i2ii22i^i^i2£i21i2ii2S the magnitude 
of the effective temperature for the given pair (T, F) does 
not indicate any dependence on c. 

Therefore one can conclude that for the GDM, i.e., for 
a disordered system with a Gaussian distribution of en- 
ergies and without any correlations between spatial po- 
sitions of hopping sites and their energies, the concept 
of the effective temperature is valid. In the next Section 
we present results for the CDM, i.e. for the correlated 
system. 



IV. EFFECTIVE TEMPERATURE FOR 

SYSTEMS WITH SPACE-ENERGY 

CORRELATIONS 

In this Section we study the effect of the space-energy 
correlations on the interplay between the temperature T 
and the electric field F with respect to the validity of the 
concept of the effective temperature. 

In Fig. [5] the numerical results are shown for the oc- 
cupation probability for each site as a function of site 
energy for three different charge concentrations at cor- 
relation length i?corr — 2. The inset shows the same 
data on a logarithmic scale. In analogy to the uncorrc- 
lated system, at low concentrations the data can be fitted 
by a Boltzmann distribution with a parameter Teg{T, F) , 
while the data at high concentrations show a Fermi-Dirac 
shape corresponding to the same parameter Tcff (T, F) as 
for low concentration. 

Fig. [5] shows the effective temperature as a function of 
the correlation length -Rcon- The effective temperature 
strongly depends on the correlation length, thus Eq. ([2]) 
is not directly applicable to a system with space-energy 
correlations. 

Before discussing the carrier mobility in a correlated 
system one should note the following. When the space- 



energy correlations are introduced as described in Sec. [Til 
the effective width a of the energy distribution of local- 
ized states decreases. Furthermore this width becomes 
smaller with increasing i?corr- Therefore one can expect 
that with increasing i?coiT the mobility would increase 
just because of diminishing the energetic disorder due to 
the effect of correlations. In order to compare the car- 
rier mobilities in systems with different correlations and 
the same energy disorder, one should therefore rescale the 
width of the energy distribution in the correlated systems 
and bring it to the value initially ascribed to the uncor- 
related systemi^ The effect of correlations on the carrier 
mobility has already been studied by computer simula- 
tions, though in the frame of the lattice model and not 
with respect to the validity of the effective temperature.— 
In order to compare our results with those of previous 
studies we also have simulated mobilities in correlated 
systems within a lattice model along with simulations in 
the random model. 

Fig. [7] shows the mobility /i as a function of the corre- 
lation length i?corr with randomly placed sites and with 
sites on a lattice for both cases: with and without rescal- 
ing (normalization) of the width of the energy distribu- 
tion. As expected, without rescaling the mobility always 
increases with the increase of the correlation length i?corr- 
However in the rescaled system there is an apparent dif- 
ference between the lattice model and the random model 
with respect to the dependence of the mobility on the 
correlation length. While in the lattice model the mo- 
bility monotonously increases with i?corr in accord with 
the results of prevous studies, ^"^ in the random model the 
dependence of /i on i?corr appears non-monotonous. At 
-Rcorr < 1 mobility decreases with i?corr- The reason for 
the low mobility at i?corr ~ 1 could be the following. The 
normalization (rescaling) procedure creates a small num- 
ber of sites with very low energies. While the general 
shape of the DOS for the correlated system is Gaussian, 




and the variance of the DOS is normahzed to be equal 
to the variance of the uncorrelated system, the DOS of 
the correlated system has longer tails, and the low-energy 
tail greatly reduces the carrier mobility. The sites in the 
tail are those sites that initially had a low energy, and 
have no neighbors within the distance Rcom- Thus they 
are unaffected by the averaging but are still rescaled with 
the normalization factor, which is typically close to 2 for 
-Rcorr = 1 ■ This idea is supported by data shown in Fig. [S] 
The dip in the mobility disappears when the charge car- 
rier concentration is increased, and the deep states are 
filled. 



V. CONCLUSIONS 



FIG. 6: The effective temperature as a function of the corre- 
lation length, for randomly placed sites, with a = 0.5. 
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FIG. 7: Mobility as a function of correlation length, with 
and without rescaling of site energies, for T = 0.4, F = 0.4, 
a = 0.5, L = 30. 



The concept of the effective temperature has been 
shown to be valid for a system with Gaussian DOS with 
and without space-energy correlations in the distribu- 
tion of localized states. From the numerical results for 
uncorrelated systems one can conclude that the effec- 
tive temperature is described by Eq. ^ with parameters 
/3 = 1.54 ± 0.2; 7 = 0.64 ± 0.2. Remarkably the valid- 
ity of the concept of the effective temperature with very 
close numerical parameters has recently been reported 
on the basis of experimental study on the inorganic dis- 
ordered material a-Si:HJ? This quantitative agreement 
between experimental data and the results of our calcu- 
lations could mean that in the material studied in"*" the 
distribution of localized states can be described by the 
Gaussian disorder model with uncorrelated distribution 
of localized states. Note that the effective temperature 
in our numerical calculations has been proven to be in- 
dependent of the charge carrier concentration. 
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FIG. 8: Mobility as a function of correlation length, with 
randomly placed sites and rescaling of the site energies, for 
different concentrations of charge carriers; T — 0.4, F — 0.4, 
a = 0.5, L = 30. 
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